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ABSTRACT 

We investigate the pc-scale kinematics and kpc-scale radio morphology of the quasar 
PG 1302-102, which may harbour a sub-pc separation supermassive binary black 
hole system at its centre as inferred from optical variability. High-resolution radio 
interferometric measurements obtained with the Very Long Baseline Array (VLBA) in 
the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) 
programme at 15 GHz at 20 epochs spanning 17 years were analysed to investigate the 
pc-scale radio structure. Archival observations with the Very Large Array (VLA) at 
1.4 GHz and 5 GHz were obtained to study the kpc-scale morphology. We find that the 
pc-scale jet is inclined within ^ 2.2° to the line of sight and has a half-opening angle of 
~ 0.2°. The parameters derived from the pc-scale radio jet are qualitatively consistent 
with those obtained from the analysis of the optical light curve of PG 1302-102. We 
obtain at least 0.08 for the mass ratio of the two black holes in the system. We find 
some indication for a helical jet structure on kpc-scale, but the directions of the inner 
and the extended radio jets are significantly different, obstructing a straightforward 
connection of the pc- and kpc-scale jets within the binary scenario. 

Key words: galaxies: active - quasars: supermassive black holes - quasars: individual: 
PG 1302-102 - radio continuum: galaxies - techniques: interferometric 


1 INTRODUCTION 

The quasar PG 1 302-102 came into focus since 
iGraham et al.l il2015l ) presented its optical variability 
with a period of 5.2 ±0.2 yr. This was found in a systematic 
search for variable sources in the V-band light curves of 
the Catalina Real-Time Transient Survey (CRTS). They 
interpreted this period as a consequence of the orbital 
motion of a supermassive binary black hole (SMBBH), and 
suggested that PG 1302-102 is a tight SMBBH candidate 
with a separation of about 0.01 pc. PG 1302-102 is a 
bright, flat-spectrum radio quasar at redshift a = 0.278 
llMarziani et al.l Il996l ) hosted by an elliptical galaxy, as 
typical for radio-loud quasars. 

The existence of SMBBHs is a natural consequence of 
the hierarchical galaxy evolution where the galaxies and 
their central black holes (BH ) grow via major mergers (e.g. 
iKauffmann fc HaehneltiboOOh . In the first phase of the evo- 
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lution, field stars carry away the energy of the binary black 
hole (BBH) and the binary shrinks under the control of 
the dynamical friction. After the depletion of stars in the 
nuclear region, an intermediate phase follows. The binary 
spends the majority of its lifetime at separations of 0.01-1 pc 
llBegelman et al.lll980h . In the third phase, the gravitational 
radiation takes over in the energy and orbital angular mo¬ 
mentum dissipation, gradually leading to a merger. The 
sub- pc separation of th e putative SMBBH in PG 1302- 
102 (iGraham et al.ll2015h suggests that this binary already 
evolved into the third, inspiral phase of the merger. For a 
recent review on the compact SMBBHs fro m observational 
point of view see lKomossa fc Zensud ll2015l) . 

Gravitational waves emitted by sub-pc separation 
SMBBHs would be observable with Pulsar Timing Array s 
(PTA; e.g. ISesana et al.ll200^ : [Manchester fc IPTAI 120131 ) . 
with sensitivity on the low-frequency regime of gravitational 
waves. Identification of possible targets is essential to con¬ 
strain the parameter field of the gravitational waveform in 
order the detect them. Based on the current parameter es- 
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timates, PG 1302-102 unfortunately remains by a factor of 
^ 10 below the detect ion threshold of even the future PTAs 
dD’Orazio et ah I [2^. 

PG 1302-102 belongs to the subclass of radio-loud 
quasars representing o nly about ten per cent of the optically 
ident ified quasars le.g. iKellermann et al.lll989l : IStocke et alJ 
I 1992 II . Studies suggest that a fast rotating BH mi ght give 
rise to relativistic jets emanat ing from the core ( Penrosa 
1 19691 : iBlandford fc Znaieklll97^ 1 , causing the radio loudness 
of the source via synchrotron radiation of the charged parti¬ 
cles spiralling in the magnetic field with relativistic speeds. 
It is plausible to suspect that BHs in radio-loud quasars 
have high spin. Probably the SMBH spin is the key factor 
to explain th e dichotomy between radio -quiet and radio-loud 
quasars fe.g. iTchekhovskov et al.1l2010l l. 

The detailed analysis of the radio jet morphology in 
active galactic nuclei (AGN) has the potential to unveil the 
properties of the central engine. The technique of Very Long 
Baseline Interferometry (VLBI) is successfully applied to in¬ 
vestigate the inner structure of the jets with the finest reso¬ 
lution. The Monitoring Of Jets in Active galactic nuclei with 
VLBA Experiment^ (MOJAVE) offers calibrated VLBI vis- 
ibilty data of jets associated with AGN at the observing fre¬ 
quency of 15 GHz. PG 1302-102 was also monitored and 
calibrat ed, and its data are available from the MOJAVE 
Survey dLister et al.ll200S^ L 

Aiming to check whether the Very Long Baseline Ar¬ 
ray (VLBA) observations reveal hngerprints of the BBH on 
the jet structure and motion on milliarcsecond (mas) angu¬ 
lar scale, we decompose the brightness distribution of the jet 
into Gaussian components at 20 different epochs over a time 
interval of more than 17 years, and investigate their long¬ 
term motion. We also use two archival Very Large Array 
(VLA) data sets of PG 1302-102 at 1.4 and 5 GHz observ¬ 
ing frequencies, to investigate the jet morphology at arcsec- 
scale. A slow precession of the spin of the BH responsible for 
launching the jet may manifest itself in a helical structure 
of the large-scale radio jet. 

In this paper, we investigate and interpret the radio 
structure of PG 1302-102 both on pc- and kpc-scales. The 
paper is organized as follows. In Section [2l we describe the 
brightness distribution model fits to archival VLBA and 
VLA observations. In Section O we discuss the morphol¬ 
ogy and kinematics of the jet of PG 1302-102. In Sec¬ 
tion 21 we interpret our results and constrain the param¬ 
eters of the SMBBH. In Section [5] we give concluding re¬ 
marks and a summary. A fiat ACDM cosmological model 
is adopted throughout the paper, with Hubble constant 
Hq = 67.8 km s~^ Mp c~^ and matter density flm = 0.308 
dPlanck GollaborationI l2015l L In this model, the angular 
scale is 4.352 kpc arcsec“^ at the redshift of the source. 


2 RADIO INTERFEROMETRIC DATA 

2.1 Gaussian decomposition of the archival VLBA 
data 

Calibrated VLBI visibilities of PG 1302-102 are available 
from the MOJAVE database llLister et al.l 1200^ 1. observed 

^ http://www.physics.purdue.edu/astro/MOJAVE/ 


in the period from 1995.57 till 2012.59 at 20 epochs. We 
used these data to reveal the radio structure at pc-scale. 
The stan dard modelfit proced ure in the Caltech Difmap 
program llSheoherd et al.lll994l l was used to decompose the 
brightness distribution of the jet into circular Gaussian com¬ 
ponents at each individual epoch. The fitted parameters of 
the circular Gaussians are the integral flux density, the po¬ 
sition of the component, and its diameter (full width at half 
maximum, FWHM). The error estimates for the fitted pa¬ 
rameters were obtai n ed in the same way as described in 
details in iKun et al.l (l2014l L Component identification was 
performed using the modelfit results, in a way that only 
slight shifts were allowed in the component parameters be¬ 
tween subsequent epochs. The brightest component at all 
epochs but 1999.55, when a newly emerged component be¬ 
came temporarily brighter, was identified with the “core” 
and was considered as the reference point. The right ascen¬ 
sion and declination of the jet components are given relative 
to the position of the VLBI “core” (Fig.[TJ. The component 
separations from the “core” as a function of time are plotted 
in Fig. 21 their relative right ascensions and declinations at 
all epochs together are shown in Fig. [3] 

2.2 Calibration and Gaussian decomposition of 
the archival VLA data 

Deep 1.4-GHz VLA observations of PG 1302-102 were ob¬ 
tained from the US National Radio Astronomy Observatory 
(NRAO) archivfl The data were taken in the most extended 
A configuration of the array, with a bandwidth of 100 MHz, 
on 1992 October 9-10 (project code: BC012). The total on- 
source integration time was nearly 1 h. Phases and ampli¬ 
tudes were calibra ted in the Ast ronomical Image Process¬ 
ing System (AIPS,[g reiserj 1^00311 in a standard way, except 
that none of the primary VLA flux density calibrators were 
scheduled in the experiment. Therefore the source OQ 208 
(J1407-I-2827) was used to define the flux density-scale, by 
adopting the total flux density value of 817 mjy obtained 
from the NRAO VLA Sky Survey (NVSS, ICondon et al.l 
Il998ll . OQ 208 is k nown to be a compact but somewhat 
variable source fe.g. IWu et al.l[2013l l. Because of the uncer¬ 
tainty in the absolute flux density calibration, we assume 
that the image brightness scale can be offset by up to ~15 
per cent. However, this does not affect our results, since we 
use the component positions in the brightness distribution 
model but not the flux densities in our subsequent analysis. 
The 1.4-GHz VLA im age of PG 1302- 102 shown in Fig. |4] 
was made in Difmap JShepherd et al.lll994l l from the cali¬ 
brated visibility data exported from AIPS, using a standard 
hybrid mapping process. 

We also obtained archival 5-GHz VLA A-configuration 
data taken on 1992 November 18 (project code: AG361). 
The total bandwidth was 100 MHz. The source PG 1302- 
102 was observed for 35 min. Phases and amplitudes were 
again calibrated in AIPS. The primary VLA flux density 
calibrator 3C 286 was used to define the amplitude scale. 
Imaging and model-fitting were performed in Difmap. The 
5-GHz image (not reproduced here) shows qualitatively the 
same two-sided arcsec-scale jet structure as the 1.4-GHz one, 

^ http://archive.nrao.edu 
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Figure 1. Temporal variation of the PG 1302-102 jet ridge line across the observing epochs of the MOJAVE radio observations with the 
VLBA at 15 GHz. Relative right ascension and relative declination are shown on the horizontal and vertical axes, respectively. The VLBI 
core lies at the (0,0) position. The black dots denote the positions of the peak of the Gaussian components at a given date indicated in 
the bottom-left corner of each panel. The light dots represent all components at all epochs in each of the panels. 
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Figure 2. Component separations from the core as a function of time. The components identified through different epochs (Cl, C2, C3, 
C4, C5, C6) are shown with different symbols. The dashed lines indicate their fitted linear proper motions. A few unidentified components 
not used for the proper motion determinations are marked with open up-pointing triangles. 




Figure 3. Relative right ascension and relative declination of the components from Cl to C6 (see in Fig.O, at all epochs. The exploded 
view of the innermost part is in the right panel. The estimated jet opening angle is marked with the dotted lines. 


© 0000 RAS, MNRAS 000, [T]-?? 





SMBH binary in PG 1302-102 5 


also known from the literature fe.g. iHutchings et al.lll9l^ . 
The extended jet components are less pronounced here than 
at 1.4 GHz, due to the higher angular resolution and the 
steep spectrum. 


3 RADIO STRUCTURE OF PG 1302-102 

While VLBI is able to map compact features and their 
variations in radio sources, its typical resolution is still far 
from being sufficient to directly image such a tightly bound 
SMBBH that is exp ected in the centre of PG 1302-102 
l|Graham et ahlEoiiTf l . Therefore we focus on an indirect way 
to identify the presence of SMBBH, and analyse long-term 
radio data to investigate the structure and kinematics of the 
jet. 


3.1 PG 1302-102 on pc scale 

Long-term VLBI monitoring data offer a way to check the 
validity of the scenario of a sub-pc separation SMBBH as 
being the central engine of PG 1302-102. Parameters of 
Gaussian components htted to archival VLBA data from 
the MOJAVE survey span 17.02 years in 20 epochs. This 
period is long enough, and the time sampling is sufficiently 
dense to allow us to investigate the jet kinematics. 

To see whether the jet ridge line is helically disturbed, 
we plot the relative positions of the Gaussian components in 
Fig.[T] If we look at the jet as a whole, its ridge line is rather 
straight and not curved. The jet points nearly to the same 
direction across the MOJAVE epochs. The average position 
angle of the jet is 0pc = 31? 6 ± 0?6 (position angles are 
conventionally measured from north through east). 

We identified 6 components of the brightness distribu¬ 
tion of the jet and followed their motions through multiple 
epochs. Their core separations versus time are plotted in 
Fig. All components show apparent superluminal motion, 
suggesting small inclination angle and large bulk speed of 
the pc-scale jet. In the linear approximation, the components 
show different proper motions. The maximum proper motion 
(i.e. the steepest slope among the lines in Fig. belongs to 
the component C4, = 0.41 ±0.04 mas yr“^. It is trans¬ 

lated to the maximum apparent speed of 3^^ = 7.5 ± 0.7, 
expressed in the units of the speed of light c. Then the min¬ 
imum bulk Lorentz factor is constrained as = 6.8, and 
the minimum jet velocity as = 0.989. 

We plot the relative right ascension and declination of 
the components collectively in Fig. [31 Interestingly, the com¬ 
ponents seem to be launched in slightly different directions, 
but after about the first 1 mas they become collimated with 
each other. The collimation is the most prominent for the 
two components at the opposite edges of the jet. Cl and 
C5. We use the outermost positions of Cl and C5 compo¬ 
nents (within ~1.5 mas separation from the core) to de¬ 
fine the opening angle of the jet cone as shown by the 
dotted lines in Fig. (3] The apparent half-opening angle is 
V-app = 5?15±0?75. 

A consequence of the Doppler boosting in a jet being 
inclined close to our line of sight is that one measures en¬ 
hanced apparent brightness temperature Tb instead of the 
intrinsic brightness temperature Tint. In case of VLBI com¬ 
ponents, the brightness temperature is calculated as (e.g. 


ICondon et al.l[l98^ : 

Tb = 1.22 X 10^2(1 ±^)^K, (1) 

where S (Jy) is the flux density, 6 (mas) is the FWHM di¬ 
ameter of the circular Gaussian component, and / (GHz) is 
the observing frequency. The apparent and intrinsic bright¬ 
ness temperatures are related to each other as Tb = J2jnt, 
where 

r(l —/Ijcost) 


is the Doppler factor, Pj is the jet speed in the units of c on 
pc scale, F = (1 — /3?)“° ® is the Lorentz factor, and 4 is the 
inclination angle with respect to the line of sight. 

We calculated the apparent core brightness tempera¬ 
tures from the fitted Gaussian component flux densities 
and diameters at each of the 20 epochs of MOJAVE data. 
Assuming the equipartition brightness temperat ure Tea ~ 
5x 10 ^° K as the intrinsic brightness temperature (iReadheadI 
Il994h . we obtained the Doppler factors. For the following, 
we adopt the median value, 5 = 18.5 ± 3.6, as the Doppler 
factor characteristic to the source. 

Then 4 and /3j can be expressed from Eq. (2] and also 
from the definition of the apparent speed 


a _ Pj sin 4 

Papp — I/O 

1 — Pj cos 4 


(3) 


Substituting the measured Doppler factor S and the proper 
motion in the jet P^^^, we obtain Pj = 0.996lg;QQ2, F = 
10.8li'g, = 2?2 ± 0?5. Then the jet inclination angle 

and the apparent half-opening angle determine the intrinsic 
half-opening angle of the jet cone as (" = V'app sin 4 = 0?20 ± 
0?05. 

The compact pc-scale radio structure of PG 1302- 
102 revealed by VLBI observations is a result of rela- 
tivistic beaming and therefore the jet inclination is small. 
[Graham et al.l (l2015tl htted the optical light curve of the 
PG 1302-102 assuming that its variability is caused by the 
periodically varying Doppler boosting in the jet. They ob¬ 
tained the following best-htting parameters: Lorentz fac¬ 
tor Fopt = 5.4 ± 0.1, intrinsic half-opening angle ^opt = 
0?5 ± 0?1, inclinatio i i ang le 4opt = 5?0 ± 0?2. Considering 
that [Graham et al.l (1201511 used optical data with an as¬ 
sumed value of the spectral index, and our results are esti¬ 
mated independently from high-resolution MOJAVE radio 
data adopting a commonly used value for the intrinsic tem¬ 
perature of the VLBI core, the two sets of parameters are 
qualitatively consistent with each other. We can thus con¬ 
clude that the pc-scale radio st ructure is compatible with 
the SMBBH model proposed bv [Graham et all ([2015[j . 

The analysis of the pc-scale morphology reveals that 
the jet length in PG 1302-102 is increasing, as seen from 
the time series plot of the jet components in Fig. [T] This 
suggests continuously strengthening jet activity. At the last 
MOJAVE epoch (2012.59), the total flux density of the jet 
components was 1.453 Jy, almost three times higher than the 
average value before this epoch, indicating a radio outburst. 


3.2 PG 1302-102 on kpc scale 

The VLA image and the fitted Gaussian brightness distri¬ 
bution model components of PG 1302-102 (Fig.[4|| are rem- 
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Figure 4. The 1.4-GHz VLA image of PG 1302-102, overplotted 
with the conical helix fitted to the components of the northern 
jet only. The dashed curve to the south of the core indicates the 
trajectory of a hypothetical jet which would be symmetric to the 
helical model fitted to the northern part. The real structure of the 
southern jet is markedly different from what would be expected 
from this simple model. The circles and the central ellipse indicate 
the Gaussian model components fitted to the visibility data in 
Difmap. Their sizes represent the FWHM diameter. The image 
peak brightness is 543 mjy beam“^. The lowest contour level is 
±4.3 mJy beam“^, further positive contours increase by a factor 
of 2. The off-source rms noise is 0.08 mjy beam“^. The restoring 
beam size (1'.'66 X is shown in the bottom-left corner. 

iniscent of a helical structure in the northern jet, and a bent 
trajectory in the southern jet. At 1.4 GHz, the “core” and 
the northern jet is represented by 7 components, while the 
southern jet is described with just two components. The lat¬ 
ter do not allow us to fit any meaningful geometric model. 
On the other hand, by assuming that the positions of the 
Gaussian components in the northern jet follow a helical ge¬ 
ometry, we could fit a si mple precession m odel to the data 
(see rotation matrices in iKun et al.l[20l3j . For this fit, the 
positions of the northern jet components from the 5-GHz 
VLA data were also used. This model provides estimates for 
the orientation and the half-opening angle of a helical jet. 

The best-fitting parameters were obtained by itera¬ 
tively performing non-linear least squares estimates using 
the Levenberg-Marquardt algorithm, such that the was 
minimized during the process. This way we got the posi¬ 
tion angle of the jet axis as 0kpc = 20? 15 ± 0?26, the in¬ 


clination of the jet axis with respect to the line of sight as 
tkpc = 65?5 ± 2?2, and the intrinsic half-opening angle of 
the jet as 4''“* = 6?5 ± 1?0. The axial growth of the jet 
during one spatial period is a = 6?'0 ± 0?H, and the radial 
growth during one spatial period is b = Cf.'T ± 0?H (reduced 
= 4.6). Our results are qualitatively consistent with the 
two-sided appearance of the kpc-scale jet which also suggests 
large inclination angle. 

We show the best-fitting northern helical jet plotted to¬ 
gether with the 1.4-GHz VLA total intensity image and the 
component locations in Fig.|4l It is important to realize that 
the components themselves are not supposed to move along 
this helix, rather their instantaneous positions form a he¬ 
lix that is expanding as the components move away from 
the core on their ballistic paths. The dashed curve in Fig. [T] 
indicates how a symmetrical helical counterjet would ap¬ 
pear in the radio image of PG 1302-102 to the south of the 
core. The structure of the southern jet is clearly different, 
indicating that this simple precessing jet model cannot be 
applied to describe it. The asymmetric kpc-scale extended 
radio structure in this object could perhaps be explained by 
jet bending due to jet-cloud interactions in the host galaxy. 


4 INTERPRETATION 

4.1 Possible precession time scale of the kpc jet 

We estimate the precession period of the northern kpc-scale 
jet, assuming that the components indeed trace out a helical 
structure (Fig.|4]). The physical expansion rate of the helical 
jet pattern is described by the deprojected axial (a) and 
radial growth (b) of the jet during one spatial period. Then 
the precession period in the rest frame of the source is 

^ VdF+lE dg (6.74 ± 1.05) X 10-* 

^ 

CPjjkpc 1 ± 2 Pjjkpc 

where /?j,kpc is the jet speed on kpc scale, and dg = 
4.352 kpc arcsec“^ is the angular scale at the redshift of 
the source. An order-of-magnitude estimate of the preces¬ 
sion period is possible if we a ssume ft.knc ~ 0.1 for the 
kpc-scale jet expansion (see e.g. IScheuenllOOSh . In this case, 
Tj is nearly 7 x 10® yr. 

4.2 The connection between the pc- and kpc-scale 
radio structures 

Contrary to the pc-scale jet, the existence of the extended 
two-sided kpc-scale radio structure in the core-dominated 
source PG 1302-102 indicates large jet inclination. How¬ 
ever, the structure of the (northern) kpc-scale jet charac¬ 
terised by our helical model (Fig. |4]) cannot be straight¬ 
forwardly explained in the context of the same binary BH 
system. The intrinsic half-opening angle of the kpc-scale jet 
(vjf'n* ~ g° 5 j is much smaller than the difference between 
the inclination of the kpc-sclale (tkpc ~ 65? 5) and pc-scale 
(t « 2?2) jets. This rules out the scenario where the pc- 
scale jet slowly processes about the symmetry axis of the 
kpc-scale jet, within its cone. In other words, assuming that 
the pc-scale jet points to the direction of the spin of the BH 
responsible for the jet launching, the kpc-scale structure can¬ 
not be interpreted simply with the precession of the black 
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hole s pin in the binary system described by I Graham et all 
(l2015h . 

However, the kpc-scale helical structure, if real, could 
possibly be reconciled with the SMB BH model if a siRnifi- 
cant jet reorientation, i.e. a spin-flip llGergelv fc BiermannI 
I 2 OO 9 I I happened in a transition phase between the produc¬ 
tion of the kpc- and pc-scale jets. We will further elaborate 
this possibility in Section 14.31 and find it unlikely based on 
the derived BBH parameters of the system. An alternative, 
albeit more speculative explanation is that the outer and 
inner structures were produced by different black holes of 
the binary system. They might have changed their activity 
status in a way that the one responsible for the kpc-scale 
jet “switched off”, and the other one now producing the pc- 
scale jet “turned on” the jet production fuelled by matter 
accretion. In this scenario, the different jet orientations are 
naturally explained by the different spin directions of the 
two central engines. 

Another, similarly speculative scenario involves the in¬ 
teraction of three BHs, where initially two BHs merged and 
reached the final coalescence. In this picture, the kpc-jet 
would show the imprint of the interaction of the two BHs, re¬ 
vealing the spin-orbit precession of the dominant spin. Then 
a spin-flip reorients the spin, causing a sudden directional 
change of the jet, and the final coalescence happens. The 
newly-merged BH catches a third BH and emanates the jet 
seen in pc scale. This jet unveils the orbital motion of the 
newly-formed binary. However, considering the estimated 
precession period of the kpc-scale jet as the consequence 
of spin-orbit precession, the gravitational lifetime w ould be 
several million years fe.g. lGergelv fc Biermannll2009l l. ruling 
out the possibility of a recent final merger in the system. 

In an alternative three-body interaction scenario, the 
third BH could have been kicked out from the system due to 
a gravitational slingshot effect, rather than being absorbed. 
This is similar to one of the prop osed explanations for the 
object CID-42 dGivoni et al.ll20ig j. In this case, the orbital 
angular momentum could have been transferred to the cur¬ 
rent primary BH, causing a spin flip without invoking a re¬ 
cent BH merger within the past several million years. Sim¬ 
ulations beyond the scope of this paper would be needed to 
prove the feasibility of this scenario. 

We may also try to explain the kpc-scale helical struc¬ 
ture assuming that it results from the precession of the ac¬ 
cretion disk of the primary BH. induced by the disk-crossin g 
motion of the secondary BH dGaproni fc AbrahamI [200^ . 
However, it would imply that the distance between the outer 
edge of the accretion disk and the primary BH is 100-1000 
times smaller than the orbital separation, ruling out the 
physical interaction between the disk and the secondary BH. 

Finally, it is also possible that the northern jet is in 
fact not helical, and therefore is not shaped by precession. 
Even in this case, the apparent discrepancy between the rel¬ 
atively large jet inclination on kpc scale on the one hand, 
and on pc scale on the other hand needs to be explained. 
Asymmetries in the dense interstellar medium surrounding 
the core, and jet-cloud interactions may r esult in the deflec¬ 
tion of the jet flow (e.g. lStocke et allll985h in a different way 
on both sides of the core. Such interactions were suspected 
by iHutchings et al.l (Il994l l who analysed the spatial corre¬ 
lation of the high-resolution optical and radio structures in 
PG 1302-102. 


4.3 Binary parameters 


We u se the BBH-in- jet model presented in our earlier 
work dKun et al.ll2014l ~) to derive parameters of the putative 
SMBBH in PG 1302-102. The underlying process revealing 
the presence of the binary is that the orbital motion of the 
jet-emitter BH modifies the average inclination angle of the 
jet launching on the orbital time scale. Using time-travel 
argument, this affects the appearance of the pc-scale jet. 
If the binary has already evolved into the inspiral phase, a 
slow spin-orbit precession may be present as well. Its period 
is much longer than the orbital period, and its signal could 
be r evealed typical l y on t he kpc-scale jet. 

iGraham et al.l (l2015l l estimated the total mass of the 
central object in PG 1302-102 from various emission lines. 
Averaging these estimates, m ~ 4 x 10® ttiq. Considering it 
as the total mass of the binary, and the res t-frame optical 
period T = 4.0 ± 0.2 yr (iGraham et al.l[2015h as the orbital 
period, the orbital separation is r « 0.01 pc. Then the post- 
Newtonian parameter e = Gmr~^c~^ « 0.002 suggests that 
the binary progressed into the inspiral phase of the merger 
{G denotes the gravitational constant). 

We can estimate the mass ratio of the two BHs, v = 
m 2 1 mi (m = mi -|-m 2 ; mi > m 2 ), using the orbital velocity 
of the dominant BH that we assume as the jet emitter: 


/ Gm'\ 2 V 
V »' / l + v' 


(5) 


The intrinsic half-opening angle of the cone along which 
the inclination angle of the jet changes is expressed as 
dKun et al.ll2014ll : 


^ = arcsin 


■i^orb cos K 
c/3j 


( 6 ) 


where k is the angle between the Newtonian angular momen¬ 
tum and the dominant spin. Substituting the expression of 
the orbital velocity from Eq. O 


ft = arccos 


/Ij sin C 1 -|- i/ 
£ 1/2 u 


(7) 


Inserting /3j ~ 0.996 and ~ 0.2° we obtained from the 
analysis of VLBA data, the mass ratio is constrained as > 
0.08. Note that if the jet inclination angle is smaller than 
2? 2, then and consequently the lower limit on the mass 
ratio would become smaller. 

Up to the second post-Newtonian order the merger dy¬ 
namics is conservative, the constants of motion being the 
total energy E and the angular momentum vector J = 
Si -I- S 2 -I- L. Here L is the orbital angular momentum. 
Si is the spin of the dominant BH, and S 2 is the spin of 
the secondary BH. The BH spin s obey precessional motion 
llBarker fc O Connelllll975l . Il979l ~) : 


Si = Qi X Si, 


(8) 


where Hi is the angular velocity of the i-th spin Si. The 
mass ratio u > 0.08 limits the ratio of the spins as 5'i52~^ < 
156 x 1 X 2 where xi and X 2 are the dimensionless spin pa¬ 
rameters of the dominant and the secondary BHs, respec¬ 
tively. In one-spin case the dominant spin Si and the New¬ 
tonian angular m omentum Ln processes arou nd J with the 
spin-orbit period (iGergelv fc Biermannll200^ : 


Tso = 


nGm - 5/2 (1 + 

c? ^ V 


(9) 
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Table 1. Binary black hole parameters derived for PG 1302-102. 


Total mass, m* {rriQ) 

ss 4 X 10“ 

Orbital period, T* (yr) 

4.0 ±0.2 

Binary separation, r* (pc) 

^ 0.01 

Post-Newtonian parameter, e 

Ri 0.002 

Mass ratio, u 

u > 0.08 

Spin-orbit precession period, Tgo (yr) 

< 14100 

Gravitational lifetime, Tgr (yr) 

< 7.2 X 10® 


* indicates pa rameters determined independently by 
iGraham et al.l ll2015l i. 


Then the dominant spin Si and the orbital angular momen¬ 
tum Ljv lie in the same plane, such that a = arccos{LN ■ J}, 
P = arccos{Si • J} and k = arccos{Si ■ Ln} {k. = a -\- /3). 
In the typical mass ratio range of the merging SMBHs 
which is 1 : 30 to 1 : 3, a spin flip is expected at the 
epoch when Si « Ln (iGergelv fc BiermannI 12009 ). Con¬ 
sidering 0.08 < V < 1/3, the ratio of the dominant spin 
and the Newtonian angular momentum is constrained as 
0.6 > Si/Ln > 0.1. In the case of PG 1302-102, this argues 
against the spin-flip as a viable explanation of the large dif¬ 
ference between the inclinations of the pc- and kpc-scale jets. 
However, a somewhat smaller / value, and consequently a 
lower limit on the mass ratio would already allow Si > Ln, 
and therefore a spin-flip to occur (for example with t = 1?2, 
C = 0?1, V > 0.04). 

In the inspiral stage, the approach and the eventual 
coalescence of the BHs is expected within the gravit ational 
time scale of the system (iGergelv fc Biermannll2003 '): 


rr _5Gm__4 (1 -I- v)^ 

GR nn Q ^ 

32c'^ V 


( 10 ) 


The mass ratio v > 0.08 limits the spin-orbit period to 
Tso < 14:100 yr, and the gravitational lifetime to Tgr < 
7.2 X 10® yr, respectively. We summarize the binary param¬ 
eters constrained by the observations of PG 1302-102 in 
Table [1] 


5 SUMMARY AND CONCLUDING REMARKS 

We analysed archival VLBA and VLA observations to re¬ 
veal the pc- and kpc-scale morphology of the radio jet of 
PG 1302-102, a quasar possibly harbouring a tight sub-pc 
separ ation supermassive binary black hole (iGraham et al.l 
I 2 OI 5 I I. The pc-scale jet components are launched in a cone 
with a half-opening angle of 0?20 ± 0?05. However, after 
leaving approximately the first mas from the “core”, they 
become collimated, resulting in a straight appearance of the 
average jet ridge line out to ~ 4 mas. 

In our interpretation, the half-opening angle of the jet 
is det ermined by the orbita l motion of the jet emitting black 
hole. IGraham et al.l ll2015ll found that the periodic opti¬ 
cal variability of PG 1302-102 is consistent with varying 
Doppler boosting, because of the small inclination and the 
high speed of the pc-scale jet. They gave a set of parameters 
describing a jet cone, similar to what we foun d by analysing 
radio interferometric data of the pc-scale jet. IGraham et aP 
ll2015|j concluded that the period most probably reveals the 
orbital period of the binary. Our results draw a physical pic¬ 
ture where the orbiting black hole perturbs the jet ejection 
at the footpoint of the jet. This gives rise to the periodi¬ 
cally changing jet direction and the periodical boosting of 


the optical flux. We constrained the binary parameters and 
using the VLBA data we determined a lower limit to the 
mass ratio of the two black holes {v > 0.08). 

By investigating the kpc-scale radio structure of 
PG 1302-102, we found that the northern part of the jet 
may be fitted with a helical geometric structure, suggesting 
a slow precession of the jet base. In this scenario, we esti¬ 
mated the precession period as Tj < 10® yr. However, if there 
is indeed a precession, the underlying physical process re¬ 
mains undetermined, since the large misalignment between 
the inclinations of the pc- and kpc-scale jets does not allow 
for a simple explanation of the latter within the framework 
of the binary model. The appearance of the kpc-scale jet 
in PG 1302-102, especially in the southern part, appears af¬ 
fected by the interaction of the jet with the dense interstellar 
medium. 

To reveal more about the radio structure of this pecu¬ 
liar source at intermediate (sub-arcsec) angular scales, in¬ 
terferometric imaging observation with e.g. the e-MERLIN 
array would be desirable. The analysis of mm-VLBI obser¬ 
vations would be helpful to check whether an exceptionally 
high flux-density outburst was indicated by the MOJAVE 
data at the last epoch (2012.59), when the monitoring was 
terminated. 
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